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Abstract 

We present a generic approach that allows us to de¬ 
velop a fully polynomial-time approximation scheme 
(FTPAS) for minimizing nonlinear functions over the 
integer points in a rational polyhedron in fixed dimen¬ 
sion. The approach combines the subdivision strategy 
of Papadimitriou and Yannakakis |22] with ideas simi¬ 
lar to those commonly used to derive real algebraic cer¬ 
tificates of positivity for polynomials. Our general ap¬ 
proach is widely applicable. We apply it, for instance, 
to the Motzkin polynomial and to indefinite quadratic 
forms x^Qx in a fixed number of variables, where Q 
has at most one positive, or at most one negative eigen¬ 
value. In dimension three, this leads to an FPTAS for 
general Q. 

1 Introduction 

Consider the problem 

(1.1) min{/(x) : X e P n Z"} 

for /: R" ^ R, P = {x e R" : Ax < b}, A e 
and b G Z"*. We use the words size and 
binary encoding length synonymously. The size of P 
is the sum of the sizes of A and b. We say that 
Problem dm can be solved in polynomial time if in 
time bounded by a polynomial in the size of its input, 
we can either determine that the problem is infeasible, 
find a feasible minimizer, or show that the problem 
is unbounded by exhibiting a feasible point x and an 
integer ray r G rec(P) := {x G R" : Ax < 0} such that 
/(x -I- Ar) — >• —oo as A — >■ 00. 

The main focus of this paper is Problem (HID with 
/(x) = x^Qx, where Q G is a symmetric matrix. 

Note that if Q is not symmetric, then we can replace it 
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hj Q' = which is symmetric and satisfies 

x^Qx = x^Q'x. The input size of Problem (HD with 
/(x) = x^Qx is the sum of the sizes of P and Q. When 
Q is positive semi-definite, /(x) is convex, whereas it is 
concave when Q is negative semi-definite. In the former 
case. Problem HD with fixed n and bounded P can be 
solved in polynomial time by |18) . whereas in the latter 
case by [5] (cf. Theorem 11.21 belowl. 

The computational complexity of mixed-integer po¬ 
lynomial optimization in fixed dimension was surveyed 
in 0 in, and they develop an FPTAS for maximiz¬ 
ing non-negative polynomials over integer and mixed- 
integer points in a polytope, respectively. They also 
use a weaker notion of approximation algorithm, that 
has also been used in m, called fully polynomial-time 
weak-approximation scheme. It allows the accuracy to 
depend on the maximum and minimum values of / on 
the feasible set. That is, let /max and /min be the max¬ 
imum and minimum values of / on the feasible region. 
Then a weak e-approximation x^ of the minimization 
problem satisfies 

[/(Xe) - /mini < e(/max “ /min)- 

When n <2 and / is a polynomial of degree d, it 
has been shown in m and [To] that Problem R.ll) can 
be solved in polynomial time for d < 3 and in |5| that it 
is NP-hard for d = 4. To the best of our knowledge, the 
complexity of Problem (HU) for fixed n > 3 and when 
/ is a polynomial of degree d = 2,3 is still unknown. 

We present an FPTAS for Problem HU) with 
/(x) = x^Qx and fixed n in two cases: For Q with 
at most one negative eigenvalue, and for Q with at 
most one positive eigenvalue. The numbers of positive, 
negative, and zero eigenvalues of Q from a triple known 
as the inertia of Q and can be computed in polynomial 
time (see [3] or [T6]l. 

By m, for fixed n, it is possible to decide in 
polynomial time if Problem HID with /(x) = x^Qx is 
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bounded or not. If it is unbounded, [9] gives a certificate 
of unboundedness, whereas if it is bounded, it returns 
an upper bound r on the size of an optimal solution that 
is polynomial in the input size. Thus we can replace P 
by P n {x G K." : ||x||oo < 2"^}. Henceforth we assume 
that P is bounded. 

We use a notion of approximation that is common in 
combinatorial optimization and is akin to the maximiza¬ 
tion version used in 017] for maximizing non-negative 
polynomials over polyhedra, except here, we extend the 
notion to allow for the objective to take negative values. 

Definition 1.1. Let Xmin be an optimal solution to 
Problem CID and let e > 0. We say that Xg is an 
e-approximate solution to Problem (ED if Xj is feasible 
and one of the following hold: 

L /(xmin) > 0 and /(x^) < (1 -k e)/(Xmin), 

2. /(Xmin) < 0 and /(x^) < Y^/(Xmin), 

3- /(Xmin) = 0 and /(x,) = /(x^m). 

We say an algorithm for Problem da is a fully 
polynomial-time approximation scheme if for any e > 0, 
in polynomial time in ^ and the size of the input, the 
algorithm correctly determines whether the problem is 
feasible, and if it is, outputs an e-approximate solution 
X£. 

In order to develop an FPTAS for classes of non¬ 
linear functions to be minimized over integer points in 
polyhedra, we propose a framework that combines the 
techniques of Papadimitriou and Yannakakis |22| with 
ideas similar to those commonly used to derive certifi¬ 
cates of positivity for polynomials over semialgebraic 
sets. Generally speaking, in the latter context one is 
given a finite number of “basic polynomials” /i,..., fm 
which are known to be positive over the integers in a 
polyhedron P. A sufficient condition to prove that an¬ 
other polynomial / is positive over P fl Z" is to find 
a decomposition of / as a sum of products of a sum 
of squares (SOS) polynomial and a basic function fi. 
A polynomial p(x) is SOS if there exist polynomials 
gi(x),... ,q^(x) such that p(x) = ^K^)- 

We would like to use a similar approach to arrive 
at an FPTAS. Again we work with classes of “basic 
functions”. Then, for a given /, we try to detect a 
decomposition of / as a finite sum of products of a 
so-called “sliceable function” and a basic function fi. 
Roughly speaking, sliceable functions — thanks to the 
result of [22] — can be approximated by subdividing 
the given polyhedron. 

For instance, the set of all convex functions pre¬ 
sented by a first order oracle that are nonnegative over 


P n Z" could serve as a class of basic functions, be¬ 
cause we can solve Problem EU for any member in the 
class in polynomial time when n is fixed. The nonneg¬ 
ativity assumption implies sign-compatibility, which is 
a necessary property of the set of basic functions that 
will be introduced in Section 0 Another example is 
the set of all concave functions presented by an evalua¬ 
tion oracle that are nonnegative over P fl Z". The same 
property holds true in this case. These two examples 
demonstrate that we consider not only polynomials ft, 
but also more general classes of basic functions. In fact, 
this is key to tackle the quadratic optimization prob¬ 
lem. For example, we can decompose the polynomial 
-\-y'^ — as the product of two non-polynomial func¬ 
tions: a basic function y/x^ + y^ ~ Nl and a sliceable 
function l^j. Our technique also applies, for 

instance, to the Motzkin polynomial, but to functions / 
that are not polynomials as well (see Theorem 12.1|1 . 

We postpone the discussion of what we mean by 
basic and sliceable functions to Section [71 As a conse¬ 
quence of our technique we easily derive the following 
result that in the optimization community was an open 
question for quite several years. 

Theorem 1.1. Let Q € be a symmetric matrix 

and let n be fixed. Then there is an FPTAS for 
Problem ED with /(x) = x^Qx in the following cases: 

(i) Q has at most one negative eigenvalue; 

(ii) Q has at most one positive eigenvalue. 

We present a proof for Theorem 11.11 in Section |3l 
Note that Theorem ll.il (pH) is equivalent to maximizing 
x^Qx with at most one negative eigenvalue. Together 
with known results on quasi-convex and quasi-concave 
optimization, Theorem 11.11 implies that for n = 3 there 
is an FPTAS for general Q. 

Corollary 1.1. Let Q G 2"^"' be a symmetric matrix 
and n = 3. Then there is an FPTAS for Problem (|l.ip 
with /(x) = k^Qk. 

In order to prove this corollary, we need the fol¬ 
lowing known results and notation, which are also used 
for our other results. Given a convex set C, a func¬ 
tion / is quasi-convex on C if the set {x G C : 
/(x) < a} is convex for all a G R. Convex functions 
are quasi-convex on any convex set in their domain. 
We say that / is quasi-concave if —/ is quasi-convex. 
A semi-algebraic set in R" is a subset of the form 
ULi HiLi {x G R" I /ij(x) 0} where fij: R" -)> 
R is a polynomial in n variables and is either < or 
= for i = 1,..., s and j = 1,..., (cf. [1]). For a poly¬ 
hedron P, denote by Pi the integer hull of P, i.e., the 
convex hull of P fl Z". 
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Theorem 1.2. ([T8], [5]) Tor A e Z""^” andh & Z™, 
let P = {x G M" : Ax < b} be a bounded polytope, and 
let n be fixed. 

(i) If C is a convex, semi-algebraic set given by 
polynomial inequalities of degree at most d and 
with integral coefficients of size at most I, and f 
is a polynomial which is quasi-convex on P, then 
inin{/(x) : X G P n C n Z"} can be solved in 
polynomial time in d, I and the size of P and the 
size of the coefficients of all involved polynomials. 

(ii) The vertices of Pj can be enumerated in polyno¬ 
mial time in the size of P. 

(Hi) If a comparison oracle for f is available, i.e., an 
oracle that given x, y G Z" decides whether /(x) < 
/(y) holds or not, and f is quasi-concave on Pi, 
then Problem HU can be solved in polynomial 
time in the size of P. 

Proof of Corollarv ll.li The inertia of Q can be com¬ 
puted in polynomial time (see [3] or m)- If Q has no 
negative eigenvalues, then / is convex because it is the 
sum of convex functions, so we can apply Theorem 11.21 
di]) with C = P. If Q has no positive eigenvalues, then / 
is concave, so we can can apply Theorem 11.21 (lull) . If Q 
has one negative eigenvalue, we can apply Theorem ll.il 
(P, whereas if there are two, we can apply Theorem ll.il 
( 0 . □ 

The remainder of this paper is outlined as follows: 
In Section [2] we explain basic and sliceable functions 
and present the main tool of this paper. In Section [3] 
we give some important properties of quadratics and 
prove Theorem ll.il 

2 Basic and Sliceable Functions 

Let P be a set of polytopes that is closed under taking 
subsets, i.e., such that for all P C P, P G P implies 
P G P. A set of basic functions fl-p is a class of functions 
with a common encoding structure for all elements and 
the following properties: 

(i) fl-p is closed under addition and multiplication by 

A > 0, 

(ii) each / G fl-p admits a FPTAS for Problem (II.ip 
over each polytope P G P, 

(hi) all functions in il-p are sign compatible, i.e., 
sign(/(x)) = sign( 5 (x)) for all f,g e flv, x e P, 
PGP. 

Examples of classes of functions H-p that satisfy 
the properties above are the set of sign-compatible 


convex function over P presented by a first order oracle 
(cf. [Ml [21]), sign-compatible concave functions over P 
presented by an evaluation oracle (cf. 0), and the set 
of polynomials that are negative on some polytope over 

P (cf. 017]). 

We next introduce a class of functions that can 
be combined with basic functions and preserve the 
property that they admit a FPTAS for Problem (II.ip . 
For this, we need the following notation. 

For k G Z*, c° G Z*, c^,...,c* G Z", and e > 0, 
define 

Pk(c°, c^,..., c\ e) := |x G K" : 

(1 -I- e)'=^-i < Ljfx.) <{1 + if kj > 1 

Pj(x) = 0 iikj =Q 

-(1 -h < Lj{x) < -(1 -f e)-''-' if kj < -1 

for j = l,...,z|. 

where Tj(x) := c-^ • x -|- c°. Let P(c°,c^,... ,c*,e) := 
Ukez''®k(c°, c^,..., c*, e). Note that for all c° G Z*, 
c^,...,c* G Z”, and e > 0 we have that Z" C 
P(c°, c^,..., c*, e). However, K" ^ c^,..., c*, e), 

because Pk(c°, c^,..., c*, e) is not full-dimensional if 
any entry of k is equal to zero. 

Definition 2.1. Let C C R”. A function s: C ^ 
M>o is called sliceable over C if there exist c° G lI, 
c^,..., c* G Z" and ( G Z>i such that for all e > 0 and 
for all k G Z*, we have that 

(2.2) s(x) < (1-F e)s(y) 

for all X, y G C n Pk(c°, c^,..., c^e/f) D Z". 

ITe define the size of a sliceable function s as the 
sum of the sizes o/c'^, c^,..., c* and the unary encoding 
size of (. 

It follows from Definition o that if a function s 
is sliceable over C with parameters c°, c^,..., c*, C, and 
s(x) = 0 for some x G C D Pk(c°, c^,..., cfe/Q D Z", 
then s(x) = 0 for all x G C'nPk(c°, c^,..., c*, e/C)nZ". 
Regardless of this propagation property, there exist 
large classes of sliceable functions. We present several 
examples of non-trivial sliceable functions below. To 
avoid the propagation of the zeros, we give choices of 
the parameters c°, c^,..., c*, ^ such that the zeros of 
the function lay on intersections of hyperplanes defined 
by Lj{x) = 0 for j G J for some J C 1,..., L 

The set of sliceable functions is not closed under 
additive inverses, but is in fact a proper semifield as the 
following lemma shows. 
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Lemma 2.1. Let s,r: C —>■ ]R>o be sliceahle functions 
and let A G K.>o- Then Xs, s + r and sr are sliceable 
with size polynomial in the size of s and r. Moreover, 
if s > 0, then ^ is sliceable with size polynomial in the 
size of s. 

Proof. Let s be sliceable with parameters c°, c^, ..., 

fs and r with parameters d°, d^, ..., d^'', fr- Set 
Cs+r = inax{Cs,Cr} and (^sr = 4max{Cs,^r}- We show 
that s + r and sr are sliceable with the same param¬ 
eters [c°; d°], c^,..., d^,..., d*'" and with Cs+r, Csr, 

respectively. Then for ( = Cs+r, Csr and 0 < e < 1, 

%i;k 2 ]([c°; d°], c\ ..., c'% d\ ..., d'L e/C) 

C Ski (cO, c\ ..., c'% e/C)ni?k=(d°, d\ ..., d'L e/C). 

Let x,y G B[ki;k 2 ]([c°;d°],c\...,c'%d\...,d'L 
e/C) nZ". By the containment, it follows that 

s(x) < (1 + e^)s(y) < (1 -h e)s(y) 

and 

r(x) <{1 + e^)r(y) < (1 -f e)r(y). 

For C = Cs-i-r-) by combining these, we see that 
s(x) -I- r(x) < (1 -I- e)(s(y) -I- r(y)). Hence, s -|- r is 
sliceable. Instead, for C = Csr, by combining these, we 
see that sr is also sliceable since 

s(x)r(x) <{1 + e^)s(y)(l e^)r(y) 

< (l + |)^s(y)r(y) 

< (1 -he)s(y)r(y). 

The function As is trivially sliceable the same 
parameters as s. Also, if s > 0, then it follows ^ is 
sliceable with the same parameters as s. Indeed, for 
X, y G i?k(c°, c^,..., cb, e/Cs) nZ” it holds that s(y) < 
(1 e)s(x). By dividing, we have ^ < (1 -h e)^- □ 


provided the degree and the dimension are fixed 
but this still implies that the best known bound for L 
is not of polynomial size in the input. Therefore there 
is no known bound on the size of C based on the input 
of the polynomial / that would admit an FPTAS for 
Problem (HID over P with this decomposition. 

Example. The polynomial M{x, y) = x'^y'^ + x^y'^ — 
3x^y^ -I- 1 is known as the Motzkin polynomial. While 
it is not SOS, it is known to be positive because 
-|- y^) ■ M{x,y) has an explicit SOS decomposition 
(see [231 P- 84]). In particular, for x,y ^0, we have 

V) = / 2 ! 2 ^ “ 2 ;) 2(1 -P xf 

[x^ + y-^) V 

-f x^{l - yf{l + yf + x^y^{x^ + y^ - 2)^) . 

Thus, for x,y 0, M{x,y) is a composition of slice- 
able functions and a convex function. That is, the pa¬ 
rameters c°, c^,..., c* are defined by cjx -I- + c° = 

Fj{x, y) for Fi{x,y),..., Fq[x, y) equal to x, y, 1 - x, 1 -f 
x,l — y,l + y respectively, and the convex function is 
(x^ +y^ — 2)^. Sign-compatibility is also given because 
all involved functions are non-negative. As we will see 
in Theorem 12.11 this implies that Problem (EH) with 
f{x,y) = M{x,y) admits an FPTAS over polytopes. □ 

The following lemma allows us to find approximate 
solutions by taking the best approximate solution over 
a decomposition of the feasible region. 

Lemma 2.2. Let C C R" such that C fl Z” ^ 0, and 
for j = 1, ..., m let Sj-. C —>■ ]R>o for j = 1,..., m 
and let gj- C —>■ R. Let 0 < e < 1 and e' = e/4. 
Let K C IP such that C fl Z” C IJkGif where 
Hk C R" are polytopes. Also, let iff > 0 such that 
< Si(x) < (1 + eO-bj for all x G Hk H Z” and 
j = 1,... ,m. 

Then, for 


Any function that is the sum of even powers of lin¬ 
ear forms is sliceable over R". Rational powers of slice- 
able functions are also sliceable. For instance, s(x) 
is sliceable for any sliceable function s. Moreover, any 
polynomial that is positive on a polyhedron P is slice- 
able on P. Handelman m showed that any polynomial 
/ that is positive over a polyhedron has a representation 
as /(x) = E|a|<L,aGZ:f Pc. nr=i(a* ’ ^ “ &*)“' for some 
L G Z+ and /3a > 0. It follows that / is sliceable over 
P. Unfortunately, the parameter f grows exponentially 
with L. The best bounds on L are given by Powers 
and Reznick [24], which are in terms of the coefficients 
of / and the smallest value that / attains on P. The 
smallest value is bounded by a polynomial in the input 


Xe := argmin 


: kG a: 

. 3=1 


where x^, is an e'-approximate solution to 

^g,(x)L^^ :xG CnSkflZ’^l , 
we have that x^ is an e-approximate solution to 


^5j (x)sj(x) : X G CnZ" 

.3=1 
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Proof. Let 


( m 'I 

Xmin := argmin^ ^gj(x)sj (x) ; x e C n Z” , 

let k S if such that Xmin S and define /(x) := 
^^1 gj(x)sj(x). Suppose first that /(x^,) > 0, which 

implies that 5 j(x^,) > 0 for all j = 1,..., m, because all 
functions in fl-p are sign compatible and Sj > 0. Then 

m 

f{4) = Y,g,{4)s,{4) 

t=i 

m 

i=i 

m 

< (1 + 

t=i 

m 

<(i+sTE min ) (^min ) 

t=l 

= (1 + e')^/(Xmin) < (1 + e)/(Xmin). 


The second inequality follows from approximate min¬ 
imality of x^,, which at the same time implies that 
5j(xmin) > 0 for all j = 1, ..., m, as otherwise /(x^,) < 
0. In particular, this implies that /(Xmin) > 0. 

Suppose now that /(x^,) < 0, which implies that 
5j (x^') < 0 for all j = 1,..., m because all functions in 
ri-p are sign compatible and Sj > 0. Then 


m m 

i=i f=i 

^ m 

- 2 + £/ y~!gj(Xmin)i^^ 
i=i 


< 


1 

(1 + e'Y 

1 

(1+102 


m 

Xniin)5j (Xniin) 

i=i 


/(Xmin) < 


1 

TT^ 


/(Xmin)- 


To derive those inequalities, we have again used the 
assumptions on and approximate minimality of 
x^,. In particular, the last inequality holds because 
/(xmin) < 0. The result follows from /(x^) < /(x^,). □ 


Let SiV) be the set of sliceable functions over all 
P G V with an evaluation oracle. We define 


C(flp) := 


m I 

gjSj I gj G fl-p, Sj G S{V), m finite > 

J 


We now show that there is an FPTAS for all 
functions in the class C{il'p). 

Theorem 2.1. Let V, Ll-p and C{Ll'p) be as above. 
Suppose that there is an evaluation oracle for f = 
G C{Ll-p) where gj G fl-p, Sj G S{'P), 
j = 1,... ,m. Then for every P G V, there exists an 
FPTAS for Problem (11.11) . where there size of the input 
is the sum of the sizes of gj, Sj for j = 1,... ,m, and 
P. 


Proof. Let / = J2jLi9j^j £ P G V, and let 

e > 0 . Feasibility of P D Z" can be determined in 
polynomial time by m- Thus we assume henceforth 
that P n Z’" ^ 0. 

Let the slicing parameters of Sj be c°, c],..., c -', Q 
for j = 1,..., TO and let Cmax = max^d }. Since P is 
bounded, there exists R G Z>i of polynomial size in 
the size of P and the sizes of Sj, j = 1,...,to, such 
that |c* • X + c°| < P for all X G P, i = 1,..., Ij, and 
j = 1,...,TO (see, for example, [25]). Set e' = f and 


N := 


log2(fl) 

Mg2(l-t-eV Cmax) 


< 


l0g2(P)(%S + 1) 


. Define 


Pk:=Pk([c?,...,c^], 

Cl, . . . , , . . . , C^, . . . , , e /((max) ■ 

In particular, IP C lJkGZ‘i+ •+!m .Bk. Then for K = 
{—N ,..., 7V}^j=i we have PflZ” C Pk- Notice 

that many choices of k G if are redundant since each 
non-empty set Pk is a cell or facet of a cell of the 
hyperplane arrangement in K." given by the hyperplanes 
c® • X + c° =0 and c® • x + c° = ±(1 + e'/Cmax)^ for 
k G {0,..., N}, i = 1,... ,lj, and j = 1,..., to. By 

pSj . we can enumerate the O (^{2N + 3)" ) 

cells of the hyperplane arrangement in polynomial time 
with n fixed. Hence, we can find a subset K' Q K with 
cardinality of polynomial size such that P C IJkGif' 

Let Xg be as in Lemma O For k G if', we 
determine a point x*^ G P D Pk D Z" using [12]. If 
no such point exists, then we remove k from if'. Set 
Fj := i:^Sj(x*^), computed using the evaluation oracle 
for Sj, for j = 1,...,TO. Note that satisfies the 
assumptions of Lemma 12.21 because Sj is sliceable and 
Cj ^ Cmax- Then, by Lemma [2?^ x^ is an e-approximate 
solution to Problem ([H]). 

We show that we can compute x^ in polynomial 
time. For each k G if', an e'-approximate solution x^, 
to 

min gj (x)P^ : x G P D Pk H Z” 

can be computed in polynomial time because Ll-p is 
closed under addition and multiplication by positive 
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scalars, and therefore J^JLi £ ^v- Finally, 

argmin{/(x^,) : k G K'} can be computed in polyno¬ 
mial time using the evaluation oracle of / because \K'\ 
is polynomial in the input size. □ 

Remark 2.1. The assumptions of Theorem \2.1\ can he 
slightly relaxed if we only consider functions of the form 
/(x) = 5 i(x)si(x), i.e., Uv = {A 51 : A > 0} for 
some gi : R.” —>■ R. and Si sliceahle. In fact, in this 
case no evaluation oracle for si is needed, because an 
e'-approximate solution to 

min {51 (x) : X G P n Rk n Z"} 

is also an e'-approximate solution to 

min { 5 i(x)Lj : x G P fl Pk H Z"} 

for any L\ >Q. Thus, L\ does not need to he computed 
explicitly. 

Remark 2.2. Theorem \2.1\ can he generalized from 
bounded polyhedra P to more general sets C, provided 
that for any bounded polyhedron P the feasibility and 
minimization problems on C H P HIP and the feasibil¬ 
ity problem on C C\1P can he solved in polynomial time 
in the size of P and the bit size reguired to represent 
C. The class of convex semi-algebraic sets in fixed di¬ 
mension is an example of such a class, because for fixed 
dimension the stated feasibility problem can be solved in 
polynomial time by m- 

Remark 2.3. The size of a sliceable function s is de¬ 
fined as the sum of the sizes of its slicing parameters, 
where unary encoding is assumed for the parameter f. 
With this definition, a polynomial bound on the cardinal¬ 
ity of K' in terms of the sizes of the .sliceable functions 
is derived in the proof of Theorem \2.1l This bound is 
not polynomial if binary encoding is assumed for the f 
parameters. However, in our main application of The¬ 
orem mi namely the proof of Theorem \l.ll those pa¬ 
rameters can be chosen to be constant. 

3 FPTAS for Quadratic Forms 

For positive definite matrices, it is well known that 
a Cholesky decomposition Q = LDLr^ can be easily 
computed, where P is a lower-triangular matrix and 
P is a diagonal matrix. Such a decomposition does 
not always exist for indefinite symmetric matrices |5]. 
However, there is a similar decomposition that serves 
our purpose, as described in the following remark. 

Remark 3.1. ([HI llbj l For any symmetric matrix Q G 
jny.n, invertible matrix S G and 

a diagonal matrix D G such that Q = SDS^ in 

polynomial time in the size of Q for fixed n. 


We thus have 

n 

x^Qx = '^di{Si ■ x)^ 

i=l 

where D — diag(d) is the diagonal matrix with Du = di, 
and Si denotes the i-ih. row of S. 

Remark 3.2. Since an FPTAS for c^/(x), c G Q>o, is 
also an FPTAS for /(x) as long as c is of polynomial 
size in the input size, S and D can be scaled to be in 
jnxn ^ /n fact, let c be the least common multiple of 
the denominators of all entries of S and D, which is of 
polynomial size in the size of Q, because S and D can 
be computed in polynomial time in the size of Q. Then 
c^ffx.) = x^(cS')(cP)(cS')^x and cS,cD G Z"^". 

Since S is non-singular, Sylvester’s law of inertia 
(see, for example [13] or [ID) states that the inertia of 
Q is the same as the inertia of D. Therefore D has 
exactly as many positive, negative and zero diagonal 
entries as Q has positive, negative and zero eigenvalues. 

The following lemmas will be useful for determin¬ 
ing the accuracy necessary for the numerical approxi¬ 
mations. 

Lemma 3.1. (see, e.G., (TU]) Fet m < n be nonneg¬ 
ative integers, G R, am 0, an ^ 0, 

and let X G be a nonzero root of the polynomial 

f{x) ■ Then\x\> v[\\n{\am\l{\am\ + \ai\) ■ 

i = m -\-1,... ,n}. 

The proof of the next lemma uses a standard tech¬ 
nique for converting an equation with rational powers 
to one with integer powers. 

Lemma 3.2. Letp,q,p',q' G Z with |p|, jgl, |p'|, |q'| < 
P, P,p' P 0, and let S = i^/p-\-q) — (v^ Fq'). If 5 0, 

then |<5| > 

Proof. We show that any solution 5 is a solution to 
a polynomial equation in terms of 6. Rewriting the 
equation for S by solving for ^/p and squaring both sides 
yields p = d'^-\-p'-\-q'^ — 2qq'-\-q'^-\-6{—2q-\-2q')-\-y/jy{26— 
2q -\- 2q'). Solving for \/p'{2S — 2q-\- 2q') and squaring 
both sides yields a degree 4 polynomial in 5 with integer 
coefficients that are polynomial in R. In particular, it 
can be shown the the absolute values of the coefficients 
are upper bounded by 48i?^. Since d is a solution to 
this polynomial equation, by Lemma l3.11 it follows that 
if d ^ 0, then |d| > □ 

With those technical tools, we can now prove The¬ 
orem o 
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Proof of Theorem \l.li Given Q £ Z”^", it is possible 
to compute S and D = diag(d) such that Q = SDS"’" 
in polynomial time by Bemark l3.ll Moreover, S and D 
can be assumed to be integral by Remark Suppose 
that Q has at most one negative eigenvalue, which 
implies that —Q has at most one positive eigenvalue. 
By reordering variables, we may assume that > 0 
for i = l,...,n — 1. Note that if d„ > 0, then 
/i(x) = is convex and / 2 (x) = x^(—(3)x is 

concave, so Problem (HID) with / = /i and / = /2 can 
be solved in polynomial time by Theorem 11.21 (jT]) and 
m respectively. Thus assume henceforth that dn < 0- 
Let 


S(x) 


n—1 


'i 




V 


dn\Sn * x| 


and 


s(x) 


n—1 




d^{Si ■ x)2 + ^-dn\Sn ■ x|. 

1=1 


Notice that /i(x) = g(x)s(x) and / 2 (x) = (—g(x))s(x). 
Let V := {P polytope : .P C P} and fl-p := {Xg : X £ 
R}. We will show that s is sliceable with size bounded 
by a polynomial in the size of Q, and that fl-p is a class 
of functions that we can minimize over the integers in 
polyhedra P £ P. Equivalently, we show that we can 
minimize g and —g over P n Z" for P £ P. The result 
then follows from Theorem o and Remark O 

First, we define the slicing parameters d> = Sj for 
j = 1,..., n, c° =0, and C = 1- Then for any k £ Z", 
e > 0, and any x, y £ Pk(c°, c^,..., c", e), we have 


5(x) = 


, y] di{Si ■ x)2 + ^-dn\Sn 

\ i=l 


< 


n—1 


\ X! + V-dn\{^ + e)Sn ■ y| 

\ i=l 


< (l + e)s(y)- 


Hence, s is sliceable with size polynomial in the size of 
Q. 

Second, we show that we can solve min{ 5 (x) : x £ 
P n Z"} exactly for any P £ P. The function g is 
convex on the halfspaces 5'n ■ x > 0 and Sn ■ x < 0, 
because it is the sum of a norm with a linear function. 
Even though g is convex, an exact first order oracle 
is not available because the function values can be 
irrational. We show the minimization problem can be 
solved by reformulating sublevel sets as convex semi- 
algebraic sets. 


We separately consider the problem on the halfs¬ 
paces S'n • X > 0 and • x < 0. The analysis for both 
cases is similar. Without loss of generality, assume that 

5'ra • X > 0. 

For /3 £ R, the inequality gfx.) < is equivalent 
to _ 

n—1 

\dn\ y] di{Si -xP < f 3 + \dn\Sn ■ X. 
i=l 

This is a second order cone constraint, which is well 
known to define a convex set [5] . It can be reformulated 
by recognizing that the left hand side is nonnegative 
and adding the condition that the right hand side is 
nonnegative too. With this, we can square both sides 
to get the equivalent relation 

n—1 

(3 3) E ^ 0’ 

(3 + • X > 0. 


In this representation, the set {x £ R" : g{x) < 
-y=^, S'n • X > 0} is a convex semi-algebraic set. It 
is thus possible to test integer feasibility of this set 
for /3 £ Q in polynomial time using Theorem 11.21 (P. 
Since P is bounded, and S and D are of polynomial 
size, it follows from [25] that we can compute R G Z, 
R > l|d||oo of polynomial size such that ||S'x||oo < R 
for all X £ P. Note that |i/— 5 (x)| < 2nP^ for all 
X £ P. We now perform binary search on /3 in the 
interval [—2nP^,2nP^] for integer feasibility of (I3.3I) 
with precision g := . Since these numbers 

are of polynomial size, the binary search can be done 
in polynomial time. We show now that the choice of 
precision yields an optimal solution by showing that 
it distinguishes function values of p—dn g{x) for x £ 

pnz’". 

Indeed, consider x, y £ P fl Z" and define <5 := 
V-dn{g{x) - g{y)). Setting p = \dn\ 1]”=/ dziSi ■ xf, 
q = -\dn\Sn -x, p' = \dn\Yi 7 P 1 diiSi ■ and q' = 
— \dn\Sn-y, we see that 6 = {^/p + q) — {y/p' + q')- Since 
p,p' > 0 and IpI, Ip'I, |g|, Iq'l < nR'^, by Lemma 15^ if 
<5^0, then |(5| > /r. Thus our choice of the precision 
for the binary search distinguishes function values of 
V~dng{x) for X £ P and hence we arrive at an exact 
algorithm that runs in polynomial time. Note also, 
that this implies that the precision p/R distinguishes 
function values of g{x) for x £ P. 

Lastly, we show that we can solve min{—( 7 (x) : x £ 
PnZ"} for any P £ P. Since g is convex, —g is concave 
and hence is also quasi-concave. By Theorem 11.21 
(jmj), the minimization problem can be solved exactly 
provided that we have a comparison oracle for —g. Our 
comparison oracle is realized by approximating —g to a 
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precision of /r/i?, which, as stated above, distinguishes 
function values of g. This precision can be obtained by 
approximating the square roots in g using, for instance 
binary search. See also |20| for approximation of square 
roots. 

Thus, we have shown that the set il-p is indeed a 
class of functions that we can minimize in polynomial 
time over polytopes P G V. □ 

Remark 3.3. In the proof of Theorem \l.ll there is an 
alternative method that avoids computing solutions to 
min{ —(7(x) : x G P H IP}, thus avoiding potential 

numerical issues when approximating square roots. In 
fact, let X G argmin{/(x) : x. is a vertex of (P fl Z”)/} 
and X G argmin{—g(x) : x G P H Z”}. Since —g is 
concave, its optimal value over P C\'IP is attained at a 
vertex of {PC\'IP)i. Therefore, because x minimizes f 
over the vertices of {P r\IP)i, we have f{5i) < /(x). 

Finally, we mention that for certain cases, we can 
obtain exact solutions to Problem (HU for /(x) = 
x^Qx where Q has at most one negative or at most 
one positive eigenvalue. For this, we use the following 
lemma. 

Lemma 3.3. Let S G Z"^" be non-singular and D = 
diag(d) with d G x Z<o. Then /(x) := x^SDS'^x 
is quasi-convex on the sets Ci := {x G K." : /(x) < 
0, 5'„ • X > 0} and C 2 '■= {x G K” : /(x) <0, Sn ■ x < 
0}. 

Proof. We show this for the set Ci ; the proof for the set 
C 2 is similar. We only need to consider the case where 
S is the identity matrix, because convexity is preserved 
under linear transform ations. _ 

Cl = {x G K" : \jYJiZl di{XiY < y/-dnXn, Xn > 
0} is a convex because it is a second-order cone [2]. In 
particular, its interior, which is given by int(Ci) = {x G 
K" : /(x) <0, Xn > 0}, is convex. 

Note that if d, = 0, then / is quasi-convex on some 
set C if and only if the restriction of / to cci = 0 is 
quasi-convex on C fl {x G K" : a;^ = 0}. Therefore, by 
reduction of variables, assume without loss of generality 
that di Q for i = 1,..., n. 

Since di ^ 0 for i = l,...,n, Ci is a pointed 
cone and therefore does not contain any lines. Let 
x,y G int(Ci) and A G (0,1). We will show that 
/(Ax-h(l-A)y) < max{/(x),/(y)}, thus proving quasi¬ 
convexity of / on int(Ci). Let 

g{X) := /(Ax-t (1 - A)y) 

= A2/(x) + (1 - A)V(y) 

n 

+ 2X{l-X)J2d^{S^■x){S,■y), 


and note that ^ is a quadratic function in A and satisfies 
5(0) = /(x) and 5(1) = /(y). If g is convex, then 
the result follows by convexity. Otherwise, g is strictly 
concave. Assume without loss of generality that /(x) > 
/(y). We will show that g is increasing on the interval 
[0, 1], thus proving the result. 

Assume for the sake of contradiction that g is not 
increasing on the interval [0, 1] and therefore has a local 
maximum on this interval, say for A. Since g is quadratic 
and strictly concave, A is also its global maximum. By 
convexity of int(C'i), Ax -|- (1 — A)y G int(C'i) and thus 
g{X) < 0. However, int(C'i) does not contain any line. 
Therefore, there is a A G M such that Ax -|- (1 — A)y ^ 
int(C'i), implying /(Ax -|- (1 — A)y) > 0 because / takes 
the value 0 on the boundary of Ci. By maximality 
of 5(A), this implies that g{X) > 0, a contradiction. 
Therefore, / is quasi-convex on int(Ci). 

Finally, if a continuous function is continuous on 
the closure of a convex set C and quasi-convex on its 
interior, it is quasi-convex on the closure of C (see, for 
example. Theorem 2.2.12 in HD- Therefore / is quasi- 
convex in Cl. □ 

Proposition 3.1. Let Q G be a symmetric ma¬ 

trix and let n be fixed. Then there is a polynomial time 
algorithm to solve Problem (HU with /(x) = x^Qx in 
the following cases: 

(i) Q has exactly one negative eigenvalue and 

/(Xmin) < 0; 

(a) Q has exactly one positive eigenvalue and 

/(Xmin) > 0. 

Furthermore, in polynomial time, we can detect if one 
of these cases occurs. 

Proof. We suppose that P fl Z” ^ 0. As mentioned 
before, we can determine the inertia of Q in polynomial 
time. 

Suppose Q has exactly one negative eigenvalue. 
Consider the sets Ci and C 2 as in Lemma [TH Since / 
is quasi-convex on each of those sets, we can minimize 
/ over Cl n P n Z” and C2 fl P fl Z" in polynomial time 
by Theorem 1 1.2 1 (111). As Ci U C2 = {x G R" : /(x) < 0}, 
we have shown how to detect if /(Xmin) < 0 and if so, 
determine an optimal solution. 

Suppose instead that Q has exactly one posi¬ 
tive eigenvalue. We consider the decomposition in 
Lemma 13.31 for —Q = S{—D)S'^, which has exactly 
one negative eigenvalue. Since we decomposed —Q in¬ 
stead of Q, the sets Ci, C2 are Ci := {x G M" : 
/(x) >0, Pn ■ X > 0} and C2 := {x G M" : /(x) > 
0, 5'n • X < 0}. Now we compute the vertices of the 
integer hulls Pf := (P ("I {x G M" : S'™ • x > 0})/ and 
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Pf := (P n {x £ R" '■ Sn ■ ^ < 0})/, which can be 
done in polynomial time by Theorem 11.21 (Ini) . A stan¬ 
dard convexity argument shows that /(xmin) < 0 if and 
only if there exists a vertex x„ of Pj or Pf such that 
/(x.) < 0. Otherwise, all vertices of P} and Pf are con¬ 
tained in Cl U ( 72 . This certifies that P 0 Z" C Ci U C 2 , 
and hence /(x^in) > 0. If this happens, / is quasi¬ 
concave on both Cl and C 2 . By Theorem 11.21 dm)) , we 
can compute an optimal solution in polynomial time. □ 
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